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ABSTRACT 

Near-Infrared (NIR) observations are presented for five Type Iln supernovae (SN f 995N, SN f 997ab, 
SN f998S, SN 1999Z, and SN f999el) that exhibit strong infrared excesses at late times (t > 100 d). 
H- and if -band emission from these objects is dominated by a continuum that rises toward longer 
wavelengths. The data are interpreted as thermal emission from dust, probably situated in a pre-existing 
circumstellar nebula. The IR luminosities implied by single temperature blackbody fits are quite large, 
> -|^q4i-42 gj.g g-i^ g^jjjj |.j^g emission evolves slowly, lasting for years after maximum light. For SN 1995N, 
the integrated energy release via IR dust emission was « 0.5 - 1 x 10^° erg. A number of dust heating 
scenarios are considered, the most likely being an infrared echo powered by X-ray and UV emissions 
from the shock interaction with a dense circumstellar medium. 

Subject headings: supernovae: individual (SN 1995N, SN 1997ab, SN 1998S, SN 1999Z, SN 1999el) — 
dust — circumstellar matter — infrared: stars 



1. INTRODUCTION 

1.1. Background 

On 3 January 1980, MerriU (1980) detected infrared 
emission from SN 1979C that was stronger toward longer 
wavelengths. Since then similar "infrared excesses" have 
been observed in a number of other supernovae including 
SN 1980K (Telesco et al. 1981), SN 1982E (Graham et 
al. 1983; Graham & Meikle 1986), SN 1982L, SN 1982R 
(Graham & Meikle 1986), SN 1985L (Ehas et al. 1986), 
SN 1993 J (Lewis et al. 1994), SN 1994Y (Garnavich, 
Noriega-Crespo & Moro-Martin 1996), SN 1998S (Gerardy 
et al. 2000; Fassia et al. 2000) and SN 1997ab (Gerardy 
et al. 2000). These infrared excesses have generally been 
interpreted as thermal emission from dust, either forming 
in the ejecta (MerriU 1980; Telesco et al. 1981; Dwek et 
al. 1983; Elias et al. 1986; Kozasa, Hasegawa, & Nomoto 
1989a, b, 1991; Lucy et al. 1991) or lying in a pre-existing 
circumstellar medium and heated by the supernova flash 
as an "infrared echo" (Bode & Evans 1980; Dwek 1983; 
Graham et al. 1983; Graham & Meikle 1986; Lewis et al. 
1994; Garnavich, Noriega-Crespo & Moro-Martin 1996). 

In the first scenario, dust forms in the cooling ejecta 
after it reaches a critical temperature. Late time ther- 
mal emission from such dust has been invoked many 
times to explain NIR excesses observed in nova outbursts 
(Gerhz 1988), where the critical formation temperature 
is typically 1000-1200 K. While such a picture is per- 
haps the simplest explanation for the dust emission (e. g., 
SN 1987A: Kozasa, Hasegawa, & Nomoto 1989a,b, 1991; 



Danziger et al. 1991; Lucy et al. 1991), it can often be ruled 
out by energy considerations. Since the dust grain cool- 
ing time is short and dust is localized in the ejecta where 
light travel times are small, the infrared luminosity must 
be limited by the energy available in the supernova at the 
time of observation. These sources include the thermal en- 
ergy stored in the dusty gas, UV/X-ray radiation from the 
supernova, and gamma radiation from radioactive heavy 
elements in the ejecta. At late times {t > 200 d) the lumi- 
nosity of observed IR excesses becomes increasingly diffi- 
cult to explain within the dust formation scenario (Dwek 
1983; Gerardy et al. 2000). 

In light of this problem, a number of authors have sug- 
gested that the observed IR excesses are due to infrared 
dust echos. In this scenario, the emitting dust is not lo- 
calized to the ejecta, but is distributed in an extended 
circumstellar nebula. As the supernova flash propagates 
outward into the nebula, the dust is heated and then ra- 
diates thermally in the infrared. Thus the late-time IR 
luminosity is powered by the maximum-light flash which, 
even after geometric dilution over a large sphere, is poten- 
tially a larger energy source than that available to dust in 
the SN ejecta. 

Such energy considerations tended to make IR echos a 
more attractive explanation for the observed infrared ex- 
cesses. For example, Dwek (1983) found that the IR emis- 
sion from the Type IIL events SN 1979C and SN 1980K 
was too bright, and that there was not enough energy in 
these objects at late times for a dust formation scenario. 
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They were, however, able to fit the observed NIR excesses 
with an echo model. Similarly, Graham & Meikle (1986) 
were able to use a IR echo model to reproduce the observed 
IR light curve of the untyped object SN 1982E. 

On the other hand, Elias et al. (1986) observed IR emis- 
sion in SN 1985L that became brighter before fading, in 
contradiction to IR echo models that predict a plateau fol- 
lowed by fading and only rise at very early times when the 
supernova is bright. Elias et al. (1986) argued that this 
was therefore a strong case for dust formation rather than 
an infrared echo. However, Emmering & Chevalier (1988) 
published models of echos in aspherical dust clouds which 
showed that IR echos from highly flattened dust distribu- 
tions could exhibit more complicated rising and falling IR 
light curves. 

Much stronger evidence for dust formation was seen in 
SN 1987A where increased infrared emission was accom- 
panied by a corresponding decrease in optical emission, 
while at the same time emission line profiles were observed 
to shift toward the blue (McCray 1993, and references 
therein). Both effects argue in favor of dust formed with 
the supernova ejecta. Similarly, in SN 1998S, the evolution 
of the hydrogen and helium line profiles also suggested dust 
formed in the ejecta as the red-shifted side of the profile 
steadily faded while the blue-shifted side remained nearly 
constant (Gerardy et al. 2000). In this case, however, the 
thermal emission in the NIR was too bright to have come 
solely from the ejecta (Fassia ct al. 2000) suggesting that 
there may be two different dust components in SN 1998S: 
newly formed dust in the ejecta and brightly emitting dust 
in the dense GSM. 

1.2. Recent Results On Type Iln Supernovae 

The NIR excesses observed in SN 1998S and SN 1997ab 
have added new dimensions to the problem (Gerardy et 
al. 2000; Fassia et al. 2000). Both objects are classified as 
Type Iln events, a sub-class of Type II supernovae that 
exhibit multi-component "Seyfert-like" hydrogen and he- 
lium emission lines (Schlcgel 1990). This class is often as- 
sociated with very slowly decaying light curves and strong 
X-ray emission (Filippenko 1997, and references therein). 
They arc thought to be interacting with a dense and per- 
haps clumpy circumstellar wind resulting from significant 
mass-loss prior to corc-coUapse (Chugai & Danziger 1994). 

The IR emission in SN 1997ab and SN 1998S was quite 
luminous (10**^"^^ erg s^^) at late epochs. For compari- 
son, the excesses seen in SN 1979C and SN 1980K were 
first observed around day 200 and decayed exponentially 
with e-folding times ~ 100 days. Yet the strong NIR ex- 
cess in SN 1997ab was observed more than three years 
after the explosion. 

Large IR hmrinositics at late epochs strained both 
the dust formation and IR echo models, especially 
for SN 1997ab. Furthermore, the NIR spectra of 
SN 1997ab and SN 1998S were nearly identical even 
though SN 1997ab was approximately three times older. 
This seemed to imply cither that these events were coin- 
cidentally observed when they would show the same spec- 
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trum, or that the evolutionary time scale of these objects 
was quite long. These issues, plus the Type Iln classifica- 
tion of SN 1997ab and SN 1998S, led Gerardy et al. (2000) 
to propose that the IR emission was powered by the cir- 
cumstellar interaction either directly through shock heat- 
ing or via UV/X-ray emission from the shocked gas. They 
were unable to distinguish between circumstellar dust and 
newly formed ejecta dust as the emission source, but Fas- 
sia et al. (2000) used K and L' photometry from an earlier 
epoch to rule out ojccta dust as the principle NIR emission 
source in SN 1998S. 

In this paper, we present further near-infrared observa- 
tions following the evolution of SN 1998S and SN 1997ab. 
We also present NIR observations of three other Type Iln 
supernovae, SN 1995N, SN 1999Z, and SN 1999cl, which 
also exhibit large late-time IR excesses. In §2, we discuss 
the observations and data reduction procedures. Near-IR 
spectral energy distributions are presented for all objects 
in §3, along with NIR spectra of SN 1995N at four epochs 
covering nearly five years of evolution. We then consider 
several emission scenarios, including interaction shock pre- 
cursor, shock heated dust, simple IR echo, and interaction 
powered echo models. We summarize our conclusions in 
§4. 

2. OBSERVATIONS 

2.1. Supernova Dates and Distances 

SN 1995N was discovered on 5 May 1995 by Pol- 
las (1995) in MCG-02-38-017. Based on optical spec- 
tra Benctti, Bouchct, & Schwarz (1995) estimated that 
it was at least 10 months old and, following Fransson et 
al. (2001), we will adopt 4 July 1994 as the date of zero 
epoch. 

SN 1998S was discovered in NGC 3877 by Z. Wan (Li, Li 
& Wan 1998; Qiu 1998) on 2 March 1998. Following Ger- 
ardy et al. (2000), we chose the date of maximum light of 
18 March 1998 (Garnavich 1998, private communication) 
as the date of zero epoch. 

The other three supernovae were discovered cither near 
or after maximum light and we adopt the discovery date 
as zero epoch for these objects. SN 1997ab was discovered 
in an objective prism image taken on 11 April 1996 for 
the Hamburg Quasar Survey, (Hagen, Engels, & Reimers 
1997). The host galaxy is the faint dwarf irregular galaxy 
HS 0948+2018. 

SN 1999Z was discovered on 8 Feb 1999 by Schwartz 
(1999) in UGC 5608, and SN 1999el was discovered in 
NGC 6951 by the BAG Supernova Survey on 20 October 
1999 (Cao et al. 1999). 

Approximate distances to the supernovae were calcu- 
lated from the recession velocities of the host galax- 
ies, using Hq — 72 km s^^ Mpc^^ (Freedman et al. 
2001). For the three nearest objects, SN 1995N, SN1998S, 
and SN 1999el, recession velocities corrected for Virgo- 
centric infall (taken from the Lyon-Meuden Extragalac- 
tic Database^) were used. The resulting distances are 26 
Mpc (1857 km s'^), 15.5 Mpc (1114 km s'^), and 24 Mpc 
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(1705 km s ^) respectively. For the more distant objects 
SN 1997ab and SN 1999Z', raw observed redshifts of 3600 
km s-i (Hagen & Reimers 1997) and 15100 km s'^ (Jha 
et al. 1999) were used yielding distances of 50 Mpc and 
210 Mpc respectively. 

2.2. Near- Infrared Photometry 

Ncar-infrarcd and .ff-band images of SN 1998S, 

SN 1997ab, SN 1999Z, and SN 1999el were obtained from 
June 1998 to November 2000 using TIFKAM on the 2.4 m 
Hiltner telescope at MDM Observatory. TIFKAM is a 
high-throughput infrared imager and spectrograph with a 
512 X 1024 ALLADIN InSb detector. Dithered sets of on 
target images were bracketed with images of nearby empty 
fields for sky subtraction. Most nights were photomet- 
ric and photometric solutions were achieved by observing 
HST/2MASS photometric standards (Persson et al. 1998) 
at a variety of airmasses throughout the night. The data 
from non-photometric nights was calibrated by matching 
the magnitudes of field stars to images from photometric 
nights. 

J-, H-, and K-hand images of SN 1995N were taken 
during the period from July 1996 to May 2001. The ma- 
jority of observations were performed on the 1.2 m tele- 
scope at F. L. Whipple Observatory (FLWO) and with 
TIFKAM on the 2.4 m Hiltner telescope at MDM Obser- 
vatory. Additional observations were obtained using the 
NASA Infrared Telescope Facility (IRTF), Apache Point 
Observatory (APO), and IRCS on the Subaru telescope. 
In order to avoid mixing data fluxed with different pho- 
tometric systems, all SN 1995N observations were fluxed 
relative to several nearby field stars. A photometric solu- 
tion for the field stars was obtained by observing Persson 
et al. (1998) standards from MDM. Thus, the SN 1995N 
data is on the same photometric system as the data from 
the other four supernovae. 

2.3. Near-Infrared Spectroscopy 

Near- infrared spectra of SN 1995N were observed from 
April 1996 to March 2001. H- and _ftr-band spectra were 
obtained on April 4 and 5 1996 using FSPEC on the 
Multiple- Mirror Telescope (MMT), and using the Infrared 
Cryogenic Spectrometer (CRSP) on the 4m Mayall tele- 
scope at Kitt Peak National Observatory (KPNO) on 16 
June 1997 and 12 February 1998. Each night of the 
MMT data was reduced separately and then combined 
into a single spectrum. An additional if -band spectrum of 
SN 1995N was observed on 5 March 2001 using the Infrared 
Camera and Spectrograph (IRCS; Kobayashi et al. 2000) 
on the Subaru Telescope. A log of all the near-infrared 
spectroscopy is presented in Table 1. 

The spectroscopic observations were taken in many 
short exposures. Between exposures the target was 
dithered across the slit to sample the chip at multiple 
locations and provide first-order background subtraction. 
Then ID spectra were extracted using standard IRAF 
tasks. Wavelength calibration was achieved using OH lines 

* http;//www.us-gemini.noao.edu/sciops/instruments/mri/NIRIIndex.! 
s NSO Kitt Peak FTS data used here were produced by NSF/NOAO 



in the night sky background for the MMT and KPNO spec- 
tra, and using arc-lamps for the Subaru data. 

A nearby F5V star from the Bright Star Catalog (Hof- 
fleit & Jaschek 1982) was used as a telluric standard 
for the MMT and KPNO data. Subaru data was cor- 
rected for telluric absorption using nearby AV and GV 
stars selected from the Gemini Spectroscopic Standard 
Star Catalogues*^. The G-dwarf was divided by a solar 
spectrum (Livingston fc Wallace 1991; Wallace, Hinkle, 
& Livingston 1993)^ to remove stellar features. This was 
used to correct for telluric absorption in the A dwarf spec- 
trum. The stellar features were fit from the corrected A- 
star spectrum and removed from the raw A-star spectrum. 
The raw A-star spectrum with the stellar features removed 
was then used to correct the target data for telluric absorp- 
tion. (See Hanson, Conti & Rieke 1996; Hanson, Rieke, & 
Luhman 1998 for a more detailed discussion of this tech- 
nique.) The instrumental response for all the runs was 
removed by matching the observed spectra of the telluric 
standards to model spectra from the stellar atmosphere 
calculations of Kurucz (1994) 

The absolute flux levels of the spectra were set by match- 
ing the average flux in each pass band to the NIR photo- 
metric data. For two of the epochs (days 1076 and 2435), 
photometric images were obtained within a few days and 
were used to set the flux. For the other two epochs, the 
NIR photometry was extrapolated using the two nearest 
photometric points for each band. In all cases, the pass 
bands of Persson et al. (1998) were assumed. We estimate 
the flux calibration is accurate to about 20%. 

There are a few issues that we have not accounted for 
which could slightly alter the results of our analysis. For 
example, the local reddening of the supernovae is not 
known for many of these objects and thus no correction 
has been applied. We have also ignored the potential 
problems of exact photometric calibration of objects which 
have spectral energy distributions (SEDs) quite different 
from those of the photometric standards; in this case very 
red continua in the H- and iT-bands and strong line emis- 
sion in J-band. In light of these issues, we have empha- 
sized the qualitative results from our analysis and treat the 
quantitative results as approximations rather than exact 
measurements. 



3. ANALYSIS 

Table 2 lists our infrared photometric results. The ob- 
served magnitudes have been transformed into fluxes us- 
ing the zero points of Bcckwith et al. (1976) and Schultz 
et al. (2001). Near-infrared spectral energy distributions 
are plotted for each of the supernovae in Figure 1. The 
SEDs have been corrected for galactic extinction using the 
measmements of Schlegel, Finkbeiner, & Davis (1998) and 
assuming an Ry = 3.1 extinction curve. Not all of the 
data for SN 1995N is plotted. Two epochs (days 1375 and 
1785) had very poor ^f-band detections and are omitted in 
Figure 1, as are all epochs with only a ii'-band detection. 
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The overall trend is a slow evolution toward redder and 
fainter near-infrared continua. In three of the objects 
(SN 1998S, SN 1999Z, and SN 1999el) we observe the onset 
of the NIR excess. For example, in SN 1998S, the /('-band 
flux remains essentially constant from day 94 to day 355 
while the J- and i7-bands fade as if the blue photospere 
fades to reveal a much more slowly evolving red continuum 
underneath. In contrast, the iT-band flux increases during 
the blue-to-red transition in SN 1999Z suggesting that the 
IR continuum might still have been rising during this pe- 
riod. SN 1995N and SN 1997ab already exhibited strong 
NIR excesses when they were first observed. During the 
observed period, SN 1997ab slowly evolved to redder col- 
ors while SN 1995N exhibited virtually no color evolution 
over nearly 1800 days. 

3.1. Single- Temperature Fits 

For data showing clear infrared excesses, we &t H — K 
color temperatures with dust emission/absorption efficien- 
cies Qx = constant, Q\ oc A^^, and Q\ cx X~^. The results 
are listed in Table 3. We note that single-temperature dust 
emission is probably not a good description for the IR ex- 
cesses, which would generally be a superposition of many 
dust temperatures. (It would be a valid description if the 
emitting dust was confined to a thin shell or very far away 
from the SN light source so that variations would be 
small.) However, since we only have NIR data it is not 
possible to constrain the long-wavelength dust emission 
spectrum. Thus, our data is really only sensitive to the 
hottest components of the dust emission. 

The best-fit curves for Qx oc are shown in Fig- 
ure 1. We note that the J-band flux always lies signifi- 
cantly above the continuum fit by the H — K color. IR 
dust emission is only a very small component of the J-band 
flux which is likely dominated by the Pa/3 emission line. As 
a result, we ignored the J-band flux in the analysis of the 
dust emission. Also, while we have plotted the Qx oc A"-'^ 
curves, the Qx = constant and Qx oc A^^ curves also fit 
the H- and ilT-band spectra of SN 1995N. Without longer 
wavelength observations we cannot constrain the dust ab- 
sorption/emission efficiency. For the rest of the analysis 
we will assume Qx oc A~^ unless otherwise noted. 

The observed H—K temperature evolution of SN 1995N, 
SN 1997ab, and SN 1998S are plotted in the upper pan- 
els of Figure 2. SN 1995N exhibits almost no tempera- 
ture evolution, cooling only from about 900 K to about 
800 K between the first and last epochs. In contrast, 
the other two supernovae show significant color evolution. 
SN 1997ab appears to cool relatively linearly from 930 K 
to 665 K over 600 days, while SN 1998S cooled quickly 
from 1140 K to 900 K over 128 days (between day 227 
and day 355), and then more slowly down to 810 K over 
the next year. Apparently the observation of nearly iden- 
tical spectra/temperatures in SN 1998S and SN 1997ab 
noted by Gerardy et al. (2000) was merely a coincidence 
as these objects are clearly evolving, albeit over a longer 
time period than the previously observed NIR excesses in 
supernovae. 

The total infrared fluxes, calculated by integrating the 



H — K temperature fits over all wavelengths, are also listed 
in Table 3 and plotted in the lower panels of Figure 2. 
These luminosities should be regarded as lower limits on 
the total infrared luminosities, as the near-IR is only sen- 
sitive to the hottest dust and consequently a significant 
flux from much colder dust could be missed in our data. 

Although SN 1995N showed little temperature evolu- 
tion, there is significant evolution of its infrared luminos- 
ity. For the first four epochs, the luminosity seems rela- 
tively constant at about 7 x 10''^ erg s~^. Then after day 
1000 the luminosity appears to decline exponentially with 
an e-folding time « 600 d. While no plateau phase was 
observed in SN 1997ab, the infrared luminosity also ap- 
pears to decay exponentially, in this case with an e-folding 
time of fti 450 d. However, the NIR behavior of SN 1998S 
was somewhat different. The inferred IR luminosity of 
SN 1998S increased ~ 50% between day 260 and day 355, 
and then faded between day 355 and day 796 at a rate 
similar to that of SN 1997ab. 

While we interpret the IR continuum as thermal emis- 
sion from dust, it is important to estimate the extent 
that free-free emission might contribute to the NIR con- 
tinuum as there are strong recombination lines observed 
in these supernovae. Fransson et al. (2001) found that 
in SN 1995N, Ha was w 8 times stronger than expected 
for standard H I recombination due to coUisional exci- 
tation and a high optical depth for Ha. However, the 
higher order Balmer lines were consistent with Case B re- 
combination. Assuming standard nebular conditions and 
Case B recombination, the H/3 emission per unit volume is 
^t^3h0 » 1.2 X lO-^^n^ erg cm'^ s'^ (Osterbrock 1989). 
Using an abundance ratio of He to H of 0.35 by number 
(Fransson et al. 2001) the free-free emission per unit vol- 
ume e// « 3.2 X 10~^^T^/^n^ erg cm~^ sec~^ (Rybicki & 
Lightman 1979), where we have assumed a velocity aver- 
aged Gaunt factor gs ~ 1-2. For T = 10^ K, this implies 
that eff~ SeHfs- Fransson et al. (2001) report an H/3 flux 
of 8.8 X 10-15 erg s^^ cm^^ from SN 1995N on day 718, 
which translates to a free- free luminosity of L// « 2 x 10^^ 
erg s~^. Compared to the IR luminosity on day 730 of 
8 x 10'*^ erg s^^, the free-free continuum is down by two 
orders of magnitude and thus does not contribute signifi- 
cantly to the observed NIR continuum. 

3.2. Comparison to X-Ray Emission 

Comparison of the observed X-ray and infrared light 

ciirves can be used to provide clues into the origin of the IR 
emission. Both SN 1995N (Fox et al. 2000) and SN 1998S 
(Pooley et al. 2001) have been detected in X-rays and the 
observed X-ray luminosities are plotted as stars in Fig- 
ure 2. 

One possible dust heating scenario is that the IR emis- 
sion is from a precursor region in the circumstellar medium 
(CSM) just ahead of the shock front that is being heated 
by X-rays from the post shock gas. In such a situation, 
where the light-travel time is small, wo might expect the 
X-ray luminosity and the IR luminosity to be well corre- 
lated. What we find, however, is that the IR evolution 
does not follow the X-ray light curve in SN 1995N, as the 
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X-ray flux shows no signs of fading nearly 200 days after 
the IR begins to fall. Furthermore, the IR flux is signifi- 
cantly greater (as much as a factor of five) than the X-rays. 
This would require that the actual emitted X-ray flux is 
much larger and that nearly all of the X-rays are being 
absorbed by the circumstellar dust and re-cmittcd in the 
IR. However, Fransson et al. (2001) found no indication of 
strong reddening in the optical spectra of SN 1995N, sug- 
gesting that the dust extinction in the CSM is low. These 
problems argue against an IR precursor interpretation. 

We note that large IR/X-ray flux ratios are seen in su- 
pernova remnants with observed ratios varying from ^ 2 
to ~ 100 (Dwek et al. 1987). In remnants the IR emis- 
sion is from hot, dusty, shocked gas where thermal emis- 
sion from coUisionally heated dust can be the dominant 
cooling mechanism. However, the IR emission from su- 
pernova remnants is typically much cooler and it seems 
unlikely that a significant amount of dust could be col- 
lisionally heated up to 800 - 1000 K without destroying 
the dust grains. Collisional heating of dust to ^ 1000 K 
would also require prohibitively large ion densities (> 10^^ 
cm~^ for Te = 10^ K; N. Evans 2002, private communi- 
cation). Furthermore, while there could be some time lag 
between changes in the X-ray flux and changes in the IR 
flux, we would not expect to see the luminosity drop in 
the IR before it does in the X-rays as seen in SN 1995N. 

3.3. Infrared Echo Models 

On the other hand, the shape of the SN 1995N IR 
light curve is similar to the luminosity evolution of sim- 
ple infrared echo models. A number of authors have pre- 
sented models describing the infrared echos of supernovae 
in spherically symmetric dust distributions (Bode & Evans 
1980; Dwek 1983; Graham et al. 1983; Graham & Mciklc 
1986). In this paper, we shall adopt the formalism of Em- 
mering & Chevalier (1988, EC88 hereafter) which adds the 
generality of non-spherical flattened dust distributions. A 
common feature of spherical models is that they predict an 
infrared light curve that rises quickly (only while the su- 
pernova is at its brightest) and then levels off in a plateau 
before flnally decaying. 

The light curve plateau results because the supernova 
flash clears out a large dust-free cavity in the center of the 
circumstellar nebula by vaporizing all the dust. At the 
edge of this cavity is the hottest and brightest dust emis- 
sion. Outside the cavity, the dust is cooler and fainter 
as the supernova flash is spread out over a larger sur- 
face. For a pulse (delta function) supernova light curve 
that explodes at time t', the observed flux at any time t 
will come from dust lying on a paraboloid described by 
t = t' + ^(1 — cos^). As time proceeds, the paraboloids 
become wider and the vertex recedes from the supernova. 
At the time = 2ry/c, the vertex of the paraboloid will 
reach the edge of the edge of the cavity. Before this time, 
the IR luminosity, dominated by emission from the cav- 
ity's edge, will be nearly constant forming the plateau in 
the IR light curve. After the paraboloid leaves the cavity, 
the emission will come from the fainter and cooler dust 
and thus the IR luminosity will fade. A plateau is also 



expected in models with more realistic SN light curves as 
long as the SN fades in significantly less time than it takes 
for the echo paraboloid to reach the cavity wall. 

The observed IR light curve for SN 1995N is at least 
qualitatively similar to the light curves of spherically sym- 
metric infrared echos. However, the evolution timescale for 
SN 1995N is much longer than for the models that were 
used to fit infrared echos to previous IR excesses in su- 
pernovae. In SN 1979C, SN 1980K and SN 1982E, the IR 
excesses were detected a couple hundred days after maxi- 
mum light, and the IR light curve was already declining. 
The time ty at which the re-emission paraboloids reached 
the cavity wall in the model fits was typically ~ 50 - 100 
days. For SN 1995N this would have to be more like 1000 
days where the observed IR light curve begins to decay 
rapidly. This time scale is set by the radius of the dust- 
free cavity, which scales as the square root of the supernova 
luminosity. Thus, the much longer time scale would imply 
either that the dust around SN 1995N is much easier to 
vaporize or that SN 1995N is much brighter (~ 100 x) than 
those other objects. A more likely solution might be that 
some event prior to core collapse cleared out the cavity. 

3.3.1. Numerical Calculations 

To test the infrared echo scenario, we have calculated the 
total IR luminosity evolution and near-IR spectra for echo 
models based on EC88. These models assume an axis- 
symmetric dust distribution where the isodensity curves 
are ellipsoids of revolution where the dust distribution has 
a radial dependence. The circumstellar dust is as- 
sumed to have an emission/absorption efficiency with a 
power-law wavelength dependence Qx oc A~". The model 
dust cloud has an inner dust-free zone, as discussed above, 
and has no sharp outer edge. The models implicitly as- 
sume that the optical depth of the dust cloud is low and 
do not account for absorption inside the circumstellar dust 
nebula, either between the SN and the heated dust or af- 
ter re-emission. Our calculations for the total infrared 
luminosity FiYi{t), and the infrared spectrum Fx{t) were 
performed by numerically integrating Equations 8 and 15 
of EC88. 

The IR dust echo models have a mimbcr of parame- 
ters that can be varied. These include the dust emis- 
sion/absorption efficiency exponent n, the size of the inner 
cavity r„ , the dust temperature at the inner surface of the 
cavity Ty, the SN light curve, the degree to which the dust 
distribution is flattened, and for non-spherical dust distri- 
butions the angle between the axis-of-symmetry and the 
line of sight. Since the observed plateau in the IR light 
curve was suggestive of a spherically symmetric dust dis- 
tribution, we used spherically symmetric echo models for 
SN 1995N. This eliminates two of the model parameters 
and also makes the numerical calculations easier. 

In principle, if the dust-free cavity were created by the 
supernova flash vaporizing the innermost dust, the cavity 
radius r^,, the peak luminosity of the supernova, and the 
dust temperature at the inner edge of the cavity T„ would 
all be related. If the dust has a vaporization temperature 
Ty, then r^ would be the radius at which the maximum 
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luminosity of the supernova would heat the dust to Ty. 
In practice, however, a number of other details about the 
dust distribution need to be known to relate these param- 
eters (e. g., grain size, composition, and absolute values of 
the absorption/emission efhcicncy and the dust density). 
Furthermore, if the cavity is the result of pre-explosion 
evolution and not vaporization due to the SN flash, then 
no relationship would exist. Therefore, we treat Ty and 
as independent parameters. 

It is convenient to parameterize the cavity size by = 
2r^/c, the epoch at which we first observe light re-emitted 
by dust on the cavity surface directly behind the supernova 
along the line of sight. This parameter sets the timescale 
for the infrared evolution. For models where the supernova 
light curve decays quickly compared to this timescale, 
is near the "kink" in the IR light curve, where the plateau 
ends and the dust emission begins to decay. For SN 1995N, 
the kink in the IR light curve occurs near t = 1000 d (see 
Fig. 2) so we consider models with ty near this epoch. 

3.3.2. Short-Flash Echo Model 

Since we have more data on SN 1995N than on the other 
four Type Iln SNe, we use this event to test the IR echo 
scenario. The solid line in the lower left panel of Figure 2 
shows the IR light curve from a model where the input SN 
light curve is a short flash compared to the echo timescale 
ty. For the model shown, the SN light curve was expo- 
nentially decaying with an e-folding time T = ty/W and 
dust emission/absorption exponent n = 1. For short-flash 
models where the supernova light curve fades much faster 
than ty, the results are quite insensitive to the exact de- 
tails of the input light curve and a short "top-hat" light 
curve would have produced essentially the same results. 

The model IR luminosity was matched to the observed 
plateau luminosity with ty scaled to 950 d to match the 
observed kink in the IR light curve. As can be seen in 
Figure 3, the resulting light curve is qualitatively simi- 
lar to the observed IR luminosity evolution of SN 1995N 
but it fades slower after the kink. However, this is not a 
strong test of the dust echo model. The "observed" IR 
luminosity is an extrapolation of the single temperature 
continuum fit to the H — K color and thus is only sensi- 
tive to the hottest dust. In contrast, the model light curve 
shows emission from a superposition of dust over a wide 
range of temperatures. As the echo ages and the emission 
comes from farther out in the nebula, the IR emission be- 
comes more and more dominated by cooler dust that is 
missed by our near-IR photometry. 

A better test is to directly compare the observed flux- 
calibrated spectra to the evolution of the model spectrum. 
The echo model calculates spectra for each epoch in di- 
mensionless wavelength units A/A,, where A^, = kTy/h. To 
make a fit for a given set of parameters, the model spec- 
trum is matched to a single epoch of the data. The model 
is shifted vertically to match the flux levels and horizon- 
tally match the shape of the observed spectrum, essentially 
fitting Ty. The same shifts are then applied to the other 
epochs and the resulting spectral evolution can be com- 
pared to the data. 



The results of such a procedure are plotted in Fig- 
ure 3. In the upper panel we plot the observed spectra 
of SN 1995N with spectra from the short-flash echo model 
for the corresponding epochs. The model luminosity and 
Ty have been flt to the second spectrum on day 1078. The 
cavity time ty for this model is 950 d and Ty w 1100 K. 
The resulting spectral evolution compares well with the 
observations, although the model luminosity fades more 
quickly. This is in contrast to the plot in the lower left 
panel of Figure 2 where the model fades more slowly than 
the "observed" IR luminosity. Again, the reason for the 
discrepancy is that the "observed" IR luminosities are only 
extrapolations based on the H — K temperature fits and 
are therefore missing a significant component of emission 
from cool dust. Since the spectral evolution provides a 
much more direct test, we use this as our primary model 
discriminant. 

While the short-flash echo model is successful enough 
to suggest that the IR echo scenario may be correct, it in- 
volves an unrealistic SN light curve. This model describes 
an echo where most of the SN energy is released quickly, 
in this case within the flrst 100 days. This timescale, while 
typical for most supernovae, is not really a good descrip- 
tion of a Type Iln supernova like SN 1995N. SN 1995N 
was observed to have a very slowly decaying light curve 
(Baird et al. 1998), a feature seen in several other Type Iln 
events (SN 1987F, SN 1978K, SN 1988Z: Filippenko 1997; 
SN 1997ab: Hagen, Engels, & Reimers 1997). Also, bright 
X-ray emission was observed in SN 1995N at very late 
times, even beyond the beginning of the observed IR de- 
cline. Both of these facts suggest that a much longer 
timescale should be used for the light source of the echo 
model. 

A longer source light curve is also easier to reconcile with 
the very large IR luminosity at late time. A simple inter- 
polation of the "observed" infrared luminosities shown in 
Figure 2 yields a total output of about 5 x 10''^ erg in be- 
tween days 730 and 2435. Integration of the model light 
curve yields a value closer to 10^° erg. Thus the energy 
released in the infrared is signiflcantly larger than a typ- 
ical ^100 day supernova flash, but is perhaps consistent 
with the long lasting interaction powered emission seen in 
other bright Type Iln supernovae. Aretxaga et al. (1999) 
found a lower limit of 2 x 10^^ erg for the emission from 
SN 1988Z, during 8.5 years of observation. Similar emis- 
sion from SN 1995N would certainly provide enough en- 
ergy to power the observed dust emission, but would again 
call for a much longer light curve for the echo model. 

3.3.3. CSM Interaction Echo Model 

While a short-flash model is relatively successful at 
matching the observed NIR spectral evolution, the in- 
put SN light curve was not a satisfactory description of 
SN 1995N and we therefore examined echo models with 
longer supernova light curves. The late-time X-ray flux 
appears to be relatively constant thereby making a long 
"top-hat" light curve an attractive possibility. For a long 
top-hat model, the IR continuum continues to rise until 
after the SN light source shuts off. If the cutoff occurs 
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prior to t-^ then the IR luminosity levels off until after ty , 
and then begins to fade. Conversely, if the cutoff epoch is 
after ty, then the IR flux continues to rise until the light 
source shuts off and then the IR flux quickly decays. 

However, what we observe in SN 1995N is not consistent 
with such a model. The NIR flux faded very slowly prior 
to the t « 1000 d kink and fairly rapidly thereafter with 
no sign of an X-ray cutoff w 200 days after the IR kink. 
Furthermore, the qualitative results from the long top- 
hat model hold true for any fairly flat input light-curve. 
Thus the rapid decline of the IR flux after the kink forces 
us to the conclusion that most of the total emission from 
SN 1995N was released well before t = 1000 d even though 
the late-time X-ray emission was still quite bright. 

This suggests a two-component light curve which we 
model as a top-hat plus constant late-time emission. The 
results of such a model are plotted in the bottom half of 
Figure 3. For the model shown, the absorption/emission 
exponent n = 2, and the echo timescalc parameter ty = 
800 d which implies a dust-free cavity radius Vy w 0.3 pc. 
The model supernova light curve consists of a bright "top- 
hat" plateau that lasts for 400 days, followed by con- 
stant luminosity emission 20 times fainter than the bright 
plateau. The flt value of Ty « 900 K for this plot. 

Figure 3 shows that this model reproduces the observed 
IR spectral evolution at least as well as the short-flash echo 
model. However, where the short-flash model was based on 
an unsatisfying description of the SN light curve, the two- 
component model is based on a more reasonable approxi- 
mation for the UV/X-ray light curve of SN 1995N. What 
this shows is that the IR excess observed in SN 1995N can 
plausibly be explained as an infrared echo. In this case 
the echo is not powered by the flash of the supernova at 
maximum light but rather by the much more slowly evolv- 
ing emission from the circumstellar interaction. We stress, 
however, that the particular model fit shown in Figure 3 is 
probably not unique. There are enough parameters that 
one could likely fit a number of echo models to the data. 

Because we did not have nearly as much data from the 
other Type Iln supernovae, we did not attempt to fit spe- 
cific echo models for them. But the infrared behavior 
of all of these objects can probably be explained by an 
interaction echo model. The very large IR luminosities 
of SN 1997ab and SN 1999Z are easiest to explain with 
such a model. In SN 1999Z, the rising IR flux seen be- 
tween day 296 and 465 would suggest that the light source 
was still emitting a significant amount of energy in this 
time period. Since we only observed a declining IR flux 
in SN 1997ab, this would suggest that we caught the IR 
emission in SN 1997ab after ty in the fading era of the 
echo. 

The behavior of SN 1998S is not suggestive of a spheri- 
cally symmetric dust echo, but is more consistent with an 
echo from a highly flattened dust distribution. Echo mod- 
els for fiattened dust nebulae with the line of sight near the 
equatorial plane do predict a bump in the IR light curve 
near ty (EC88). This would be in line with other evidence 
for deviations from spherical symmetry in SN 1998S. Spec- 
tropolarimetry indicated that both the SN ejecta and the 



circumstellar environment were aspherical (Leonard et al. 
2000; Wang et al. 2001). Also a well-resolved triple-peak 
line profile observed for the hydrogen and helium lines of 
SN 1998S was interpreted by Gerardy et al. (2000) as emis- 
sion from interaction with dense, clumpy circumstellar gas 
in a ring or disk surrounding the supernova and probably 
observed from a fairly high inclination. 

4. CONCLUSIONS 

We have presented near-infrared photometry and spec- 
tra from five recent Type Iln supernovae that exhibited 
large late-time infrared excesses. The infrared excesses in 
these objects seem to have a much slower evolution than 
those that have been previously presented in the literature, 
remaining bright many years after maximum light. The re- 
sulting energy release via IR dust emission is quite large, 
approximately 0.5-1 x 10^° erg from SN 1995N. The long 
IR light curve and large luminosity cannot be reconciled 
with emission from newly formed dust in the ejecta. In ad- 
dition, for SN 1995N and SN 1998S, the IR luminosity is 
significantly larger than the late-time X-ray emission, rul- 
ing out emission from a shock precursor. IR emission from 
coUisionally excited dust in shocked gas also seems unlikely 
as the observed temperatures are fairly high (~ 1000 K). 
Moreover, in SN 1995N the IR fiux declines rapidly after 
day 1000, while the X-ray flux shows little sign of fading 
even « 200 days later. 

On the other hand, the data do appear to be consistent 
with an infrared echo of the bright, slowly- fading emission 
from the circumstellar interaction. Simple IR echo models 
are able to reproduce the observed near-infrared spectral 
evolution of SN 1995N and such a scenario helps explain 
the slow evolution and large IR luminosities seen in these 
objects. In principle, detailed studies of the IR cchos of 
these objects could provide important information about 
the circumstellar environment of Type Iln supernovae, and 
might lend some insight into their progenitors. In practice, 
however, to provide significant constraints the supernova 
will have to be very well observed, including early time 
UV and X-ray observations, as well as good near- and 
mid-infrared observations of the echo. 
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Table 1 



Near-Infrared Spectroscopy of SN 1995N 



Date 


JD 


Epoch'^ 
(days) 


Wavelength 
Coverage (/ini) 


Exposure 
(sec) 


Observatory/Instrument 


04 Apr 1996 


2450179 


640 


1.45-1.75 


1800 


MMT + FSPEC 








1.98-2.43 


3360 




05 Apr 1996 


2450180 


641 


1.45-1.75 


3840 


MMT + FSPEC 








1.98-2.43 


1920 




16 Jun 1997 


2450617 


1078 


1.45 1.73 


3600 


KPNO 4ni + CRSP 








1.96 2.49 


1200 




12 Feb 1998 


2450858 


1319 


1.44 1.72 


900 


KPNO 4ni + CRSP 








1.92-2.50 


900 




05 Mar 2001 


2451974 


2435 


1.93-2.49 


3600 


Subaru + IRCS 



'^Adopted date for zero epoch: 04 Jul 1994 
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Table 2 



Near-Infrared Photometry 



uate 


JL) 


rLpocn 


J 


rr 
H 


A 


(J bservatory 


CIV 1 ooE^ivr 
SIM lyyfjiN 














04 Jul 1996 


2450269 


730 


16 90 ± 03 


15 55 ± 03 


13.41 ± 0.04 


IRTF 


25 Mar 1997 


2450533 


994 


17.51 ±0.11 


15.90 ±0.05 


13.55 ±0.05 


FLWO 


25 Apr 1997 


2450564 


1025 


17.35 ±0.13 


16.07 ±0.05 


13.61 ±0.07 


FLWO 


12 Jun 1997 


2450612 


1073 


17.57 ±0.12 


15.99 ±0.07 




APO 


15 Jun 1997 


2450615 


1076 






13.64 ±0.05 


FLWO 


iy Mar iyyo 




ioOo 




iO.Oo ± U.U i 


14. zi zb U.U4 


r LiVVU 




z^touy 1^ 


1 'Xl^ 
lO / 


^ io.oi 


1 7 on _l_ n 1 
1 / .zy Hz u. ±o 


1 4 9fi -1- OQ 
i^.zD ZIZ u.uy 


FT wn 


1 A Tnn 1 QQR 


z^uuy 1 y 






1 7n -u 10 


1 A on _|_ A fizL 


WO 


Of^ TTaK 1 QQQ 
uu JrcL) lyyy 


94f^1 91 


1 R77 




1 1 .U(J ZIZ u. ± ± 


1 4 Qo _i_ n 07 

±^.OZ ZIZ U.U / 


FT WO 


9Q A/Tnv 1QQQ 
zy ividy lyyy 










1 4 QQ -1- OR 
I'rt. yy ZIZ u.uo 


FT WO 

JT J_jVV w 


11 Anr 9000 


94^1 f{Af{ 


91 07 


1 Q ^7 _|_ n 9*? 




1 c f?Q _|_ n OR 
lo.uy ZIZ u.uo 


ATT) AT 


12 Apr 2000 


2451647 


2108 




18.15 ±0.12 




MDM 


14 Apr 2000 


2451649 


2110 






15.85 ±0.10 


FLWO 


05 Mar 2001 


2451974 


2435 




18.82 ±0.06 


16.31 ±0.06 


Subaru 


30 Mar 2001 


2452000 


2461 






16.41 ± 0.05 


MDM 


01 May 2001 


2452032 


2493 






16.47 ± 0.14 


MDM 


SN 1997ab 














09 Mar 1999 


2451248 


1062 


17.58 ±0.05 


15.67 ±0.03 


13.86 ±0.04 


MDM 


18 May 2000 


2451684 


1498 




17.96 ± 0.13 


15.54 ± 0.04 


TV J'T^Tl T 

MDM 


26 Nov 2000 


2451876 


1690 




19.30 ±0.18 


16.48 ± 0.06 


MDM 


SN 1998S 














20 Jun 1998 


2450986 


94 


14.37 ±0.03 


14.07 ±0.05 


13.46 ±0.07 


MDM 


31 Oct 1998 


2451119 


227 


16.29 ±0.05 


14.83 ±0.04 


13.50 ±0.04 


MDM 


03 Dec 1998 


2451152 


260 


16.74 ±0.08 


15.05 ±0.06 


13.62 ± 0.07 


MDM 


08 Mar 1999 


2451247 


355 


17.26 ±0.07 


15.46 ±0.04 


13.57 ±0.02 


MDM 


09 Mar 1999 


2451248 


356 


17.28 ±0.08 


15.46 ±0.03 


13.59 ±0.03 


MDM 


22 May 2000 


2451688 


796 




17.57 ±0.19 


15.38 ±0.09 


MDM 


SN 1999Z 














01 Dec 1999 


2451515 


296 


17.71 ±0.05 




16.49 ±0.05 


MDM 


18 May 2000 


2451684 


465 


18.25 ±0.08 




16.07 ±0.07 


MDM 


23 Nov 2000 


2451873 


654 


> 20.11 


19.55 ±0.23 


16.84 ±0.06 


MDM 


SN 1999el 














18 May 2000 


2451684 


198 


18.33 ±0.15 




16.29 ±0.08 


MDM 


23 Nov 2000 


2451873 


387 


> 18.99 


> 18.53 


16.84 ±0.12 


MDM 



Adopted dates for zero epoch — SN 1995N: 04 Jul 1994; SN 1997ab: 11 Apr 1996; SN 1998S: 
18 Mar 1998; SN 1999Z: 08 Feb 1999; SN 1999el: 02 Nov 1999. See text for details. 
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Table 3 
Single Temperature Fits 



Epoch 




Qx oc 1 


( 


3a oc 


( 


3a oc A-^ 


(days) 


Td (K) 


Ld (10*" erg s-^)" 


Td (K) 


Ld (10*" ergs-^)" 


Td (K) 


Ld (10*" erg s-^r 


SN 1995N 














730 
994 
1025 

1073,1076 

1353 
1440 
1677 

2107,2108 
2435 


930 ± 20 
860 ± 25 
830 ± 30 
860 ± 30 
840 ± 25 
845 ± 40 
900 ± 50 
830 ± 45 
815 ± 25 


89 ± 8 
101 ± 12 
108 ± 19 
92 ± 13 
59 ± 8 
53 ± 10 
27 ± 6 
16 ±4 
10 ± 2 


830 ± 15 
770 ± 20 
745 ± 25 
775 ± 20 
755 ± 20 
760 ± 30 
805 ± 40 
745 ± 40 
735 ± 20 


65 ± 5 
71 ± 8 
75 ± 11 
65 ±8 
41 ± 5 
37 ±6 
19 ±4 
11 ±3 
7±1 


745 ± 15 
700 ± 15 
680 ± 20 
700 ± 20 
685 ± 20 
690 ± 25 
725 ± 35 
680 ± 30 
670 ± 15 


52 ± 4 
55 ± 5 
57 ±8 
50 ± 6 
32 ± 4 
29 ± 4 
15 ± 3 
8 ± 2 
5 ±0.6 


SN 1997ab 














1062 
1498 
1690 


1060 ± 20 
835 ± 50 
730 ± 55 


156 ± 11 
65 ± 14 
48 ± 17 


930 ± 15 
750 ± 40 
665 ± 45 


120 ± 7 
45 ± 9 
31 ± 10 


825 ± 15 
685 ± 30 
610 ± 40 


100 ± 5 
35 ± 6 
22 ± 6 


SN 1998S 














227 
260 
355 
796 


1350 ± 45 
1275 ± 70 
1025 ± 15 
910 ± 100 


15 ± 0.8 
14 ± 1.5 
21 ± 0.9 
5.5 ± 1.9 


1140 ± 35 
1085 ± 52 
905 ± 10 
810 ± 75 


12 ± 0.5 
11 ± 1 
16 ± 0.6 
4.0 ± 1 


990 ± 25 
950 ± 40 
805 ± 10 
730 ± 60 


11 ± 0.4 
10 ± 0.8 
13 ± 0.4 
3.1 ± 0.8 


SN 1999Z 














654 


760i«° 








630i|| 




SN 1999el 














387 
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Fig. 1.— Spectral energy distributions for SN 1995N, SN 1997ab, SN 1998S, SN 1999Z, and SN 1999el. The curves are the best fit single 
temperature blackbody curves for the observed H — K colors, assuming a absorption/emission efficiency. 
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Fig. 2. — The upper panels show the observed H — K temperature evolution for SN 1995N, SN 1997ab, and SN 1998S, assuming a 
absorption/emission efficiency for the dust. The lower panels plot the inferred IR luminosities from the H — K temperature fits (triangles) 
and observed X-ray luminosities (crosses). The solid line in the lower SN 1995N plot is the IR luminosity from an IR echo model with a short 
exponential flash as the input supernova light curve. 
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Fig. 3. — Observed near-infrared spectra of SN 1995N compared to calculated spectra from infrared echo models. For the model in the top 
panel the input SN light curve was a short exponential decay. The model in the bottom panel results from a two-level "top-hat" model with 
a long bright plateau followed by a fainter constant source. The epochs of the observed data are 640, 1078, 1319, and 2435 days, from top to 
bottom. 



